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The transport interaction rates of elastic scattering processes of thermal partons in the quark- 
gluon plasma are calculated beyond the leading logarithm approximation using the effective per- 
turbation theory for QCD at finite temperatures developed by Braaten and Pisarski. The results 
for the ordinary and transport interaction rates obtained from the effective perturbation theory are 
compared to perturbative approximations based on an infrared cut-off by the Debye screening mass. 
The relevance of those interaction rates for a quark-gluon plasma possibly formed in ultrarelativistic 
heavy ion collisions are discussed. 

PACS numbers: 12.38.Mh, 12.38.Bx 



I. INTRODUCTION 



Interaction or damping rates of parton scattering processes in a thermalized quark-gluon plasma (QGP) are of great 
physical significance. For example, the inverse rates related to the elastic scattering of thermal partons {gg — > gg, 
gq — > gq, qq — s- qq) give the mean free paths and typical interaction times i.e., relaxation times used in the collision 
term of the Boltzmann equation. Hence they can be used for an estimate of the thermalization time of a pre- 
equilibrium parton gas in ultrarelativistic heavy ion collisions Q and the maintenance of the thermal equilibrium 
^Isq . by comparing the interaction rates versus the cooling rate. Of course thermalization times calculated in this way at 
QQ ' finite temperature should be valid only for situations which do not start too far from equilibrium. On the other hand, 
this thermalization time does not depend on the somewhat ambiguous definition of the beginning and the end of the 
pre-equilibrium stage as the one deduced from numerical studies of ultrarelativistic heavy ion collisions (HIJING , 
[ parton cascade ||^). After all, both methods lead to similar results, namely a fast thermalization of the gluons i.e., an 
(~| 1 isotropic and exponential momentum distribution, after about 0.2 fm/c for LHC and 0.3 fm/c for RHIC Q. The full 
rS I equilibration of the phase space density can be investigated, on the other hand, using the inelastic interaction rates 
JL , e.g., gg —> ggg and gg —> qq, describing the chemical equilibration process of the QGP Q|. 

Furthermore the elastic rates are the basic inputs for the energy loss of a parton in the QGP [^-p^ and the viscosity 
(~| , of the QGP |]l3| . The first quantity is related to a possible signature of the QGP in ultrarelativistic heavy ion collisions 
!" ' (jet quenching) ||7|,|l4|Jl5| , while the latter provides informations about the role of dissipation in the hydrodynamic 
expansion phase [|16{. 



uantum field theory at finite temperature 
17|. In particular, the puzzle of the gauge 



Finally, damping rates are especially suited for studying problems of 
^ I like gauge dependences and infrared singularities in perturbation theory 

dependence of the gluon damping rate at rest (plasmon damping) in naive perturbation theory started a recent 
developement in finite temperature field theory, which lead to a powerful effective perturbation theory . Assuming 
the weak coupling limit {g <C 1) effective propagators and vertices are obtained by resumming self energies and vertices 
in the high temperature limit (hard thermal loops), which are used for soft external momenta of the order gT, while 
bare Green's functions are sufficient for hard momenta of the order T. In this way gauge independent results for 
physical quantities are found, which are complete to leading order in the coupling constant. (The gauge independence 
of this method is still under discussion by evaluating the damping rates in a general covariant gauge. However, using 
an appropriate infrared regularization for the gauge dependent part, independence of the gauge fixing parameter can 
be shown |[l9|-[2^.) In addition, screening effects are included yielding an improved infrared behavior. In summary, 
the Braaten-Pisarski method provides a consistent treatment for quantities which are sensitive to the momentum scale 
gT, meaning a crucial improvement compared to the naive perturbation theory at finite temperature. 

Here, we will present calculations of the interaction rates based on the Braaten-Pisarski technique and study their 
physical relevance for the QGP. Two different kinds of interaction rates are considered. In the next section we will 
discuss the ordinary interaction rate F - in the following simply called interaction rate - of a parton with hard 
momentum, defined by F = n cr, where n is the particle density of the QGP and a the cross section of the scattering 
process under consideration. In naive perturbation theory without resummation F ~ T is found (a^ — g'^/Aw). 
However, F turns out to be quadratically infrared divergent. Using the Braaten-Pisarski method, F ~ as T ln(l/as) 
follows, which is only logarithmically infrared divergent [RJ23-|2^. The interaction rate is enhanced by a factor l/us 
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compared to the one obtained from naive perturbation theory due to the infrared regularization inherent in the 
resummation method. The logarithmic term, reflecting the logarithmic infrared divergence, arises from the sensitivity 
of the rate to the momentum scale g^T^ which cannot be treated within the Braaten-Pisarski method. 

In section 3, we will discuss the transport interaction rate defined by T trans = n crtrans, where the transport cross 
section Utrans — J d,<T (1 — cos 9) enters the collision term of the transport equation in the case of a plasma with long 
range interactions, for which the interaction rate is dominated by distant collisions js^. Here 9 denotes the scattering 
angle in the center of mass system. Hence the mean free path as well as the thermalization time may rather be given 
by the inverse of the transport interaction rate |^0|. The relevance of this rate for the QGP has been pointed out in 
Ref. The transport interaction rate is also closely connected to the shear viscosity beyond the relaxation 

time approximation |p^ , ^ |. 

The factor 1 — cos 9 cuts off small scattering angles and therefore small momentum transfers, too. Thus the 
transport interaction rate is only logarithmically infrared divergent in naive perturbation theory (Coulomb logarithm 
pO|), but finite using the Braaten-Pisarski method since it is sensitive to the scale gT only. Hence dynamical 
screening of the magnetic interaction, which is included in the effective gluon propagator, is sufhcient [Q leading to 
^trans ^ o^T ln(l/as) [Il3| , p9|j3ll . Note the a1 dependence as expected from naive perturbation theory since there 
is no enhancement by a factor l/a^ caused by a quadratic infrared divergence. Also the logarthmic term stems from 
the sensitivity to the scale gT instead of g^T as it is the case for P pot] . 

Up to now the transport interaction rates based on the Braaten-Pisarski method have been calculated only within 
the leading logarithm approximation [p^ i.e., the coefficient of l/cts under the logarithm has not been determined. 
However, if we want to extrapolate the result for the transport interaction rate in the weak coupling limit to realistic 
values of = 0.2 - 0.5, we should not neglect this coefficient. 

One aim of the present paper is the calculation of the transport interaction rate beyond the leading logarithm 
approximation using the Braaten-Pisarski method. Furthermore we will discuss the consequences of P and Ttrans 
for realistic values of as. Finally, we will compare the results of the Braaten-Pisarski method with the widely used 
approximation of the naive perturbation theory regularizing the infrared singularities simply by using the Debye 
screening mass as an infrared cut-off. 



II. INTERACTION RATES 

There are two equivalent ways of computing the interaction rates either using matrix elements or self energies . 
Considering for example elastic quark-quark scattering we may find the corresponding interaction rate Tqq from the 
matrix element via ||^ 

^^'^YpJ (2^ - "^(P J J2nf2k 

r rPk' 

X J [1 - np{k')] {2iTf 5\P + K-P'-K') 6Nf {\M{qq - qq)\'). 

(1) 

The four momentum of the incoming particle is denoted by P = (p, p), where quark masses are neglected and p = |p|. 
(We consider only the case of a QGP containing up and down quarks for which the bare masses are negligible compared 
to the temperature of the QGP.) The quark of the heat bath from which the quark under consideration is scattered off 
has the momentum K = {k, k). Momenta with a prime belong to the outgoing particles. The Fermi-Dirac distribution 
functions are given by np(k) — l/[exp(fc/r) + 1], while Nf denotes the number of thermalized flavors in the QGP. 
The matrix element is averaged over the spin and color degrees of freedom of the incoming particles and summed over 
the ones of the final state. The factor 6iV/ comes from summing over the possible spin, color, and flavor states of the 
quark with momentum K. The diagrams which enter the matrix element to lowest order are shown in Fig. la. 

In the case of quark-antiquark scattering we have to replace the diagrams of Fig. la by the ones of Fig. lb. In 
the case of quark-gluon scattering we have to substitute the Fermi distributions by the Bose-Einstein distributions 
nsik) = l/[exp(fc/T) — 1], the Pauli blocking factor 1 — np{k') by the Bose enhancement factor 1 -I- nB{k'), and the 
factor 6Nf by 16 (number of gluonic degrees of freedom). If we wish to consider the scattering of an incoming gluon 
by the partons of the QGP we have to deal with the diagrams of Fig.lc and Id. 

The second possibility to determine the quark interaction rate is given by the imaginary part of the quark self 
energy on mass shell 

r,(p) = [1 - "f(p)] tr [^P, ImT.{p, p)] . (2) 
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The equivalence of the expressions (|^) and (||) can be seen from cutting the self energy of Fig. 2a through the fermion 
lines. Eq.(|^) is the starting point for applying the Braaten-Pisarski method by considering the quark self energy of 
Fig. 2b, where the effective gluon propagator contains the resummed one-loop gluon self energy in the high temperature 
limit 1^,^. (Note that Fig. 2b contains the diagram of Fig. 2a if the high temperature limit is used for the fermion 
loop, called hard thermal loop approximation |18|, in the latter.) It is not necessary to take an effective quark-gluon 
vertex into account because the external quark with a momentum of the order of the temperature is hard. Also 
because the interaction rate falls olf rapidly for large momentum transfers g = |p — p'|, F ^ / dq/q"^, i.e., only small 
momentum transfers contribute, a bare quark propagator is sufficient. An effective quark propagator and quark-gluon 
vertex would contribute to higher order in only. 

Since the final result for observables using the Braaten-Pisarski method is gauge independent we are free of choosing 
any gauge. Using Coulomb gauge and the approximation p, A: ~ 3T ^ q ^ gT the quark interaction can be written 
as i 



/ dqq — 

271- Jo J-q ^ 



pi{uj,q) + ( 1 - ) Pt{uJ,q) 



(3) 



where Cf = 4/3 is the Casimir invariant of the fundamental representation. The interaction rate of a hard gluon 
is simply obtained by replacing Cp by the Casimir invariant of the adjoint representation Ca = 3 p6[ |. The four 
momentum of the exchanged gluon is denoted by P — P' — K' — K = Q — {uj,q), and the discontinuous parts of 
the longitudinal and transverse spectral functions, related to the effective gluon propagator A/^t through p/.f(a;, g) = 
ImAi/{Lu, q)/Tr, are given by {-q < uj < q) @: 

pi{uj,q)^ 




(4) 



where = (1 + Nf/Q) may be interpreted as an effective gluon mass generated by the interaction with the 

thermal ensemble of the QGP. This expression was derived assuming the high temperature approximation w, g <C T. 

The integration over q can be performed analytically, whereas the one over uj has to be done numerically. The 
result for the longitudinal part of the interaction rate corresponding to the exchange of a longitudinal gluon is given 

by 



= 1.098 Ci. as T. 



(5) 



Note that the interaction rate is independent of the number of flavors Nf due to a cancellation of the factors 
in (Q) after integration. Although a larger number of thermal flavors corresponds to an increase of the number of 
scattering partners enlarging the interaction rate, the screening mass is also increased cancelling this enlargement. 
Also the interaction rate does not depend on the external momentum p in the p ^ uj,q limit. 

The transverse part of the interaction rate F* , on the other hand, is still infrared divergent, although the infrared 
behavior has been improved from a quadratic singularity in the bare two loop case (Fig. 2a) to a logarithmic due to 
dynamical screening of magnetic fields. On the other hand, there is no static magnetic screening in the transverse 
part of the spectral function (or the effective propagator) i.e., the denominator of pt vanishes in the static limit 
Lu ^ 0, q ^ 0, while the denominator of pi is given by = 3™^ in this limit. In other words, the high temperature 
limit of the perturbatively calculated gluon self energy contains static electric screening (Debye mass p^) but no 
static magnetic screening. In QCD, however, static magnetic screening may arise non-perturbatively from monopole 
configurations due to the gluon self interaction. Indeed, lattice as well as semiclassical calculations show the excistence 
of a magnetic screening mass, m^^g ~ 15 Thus static magnetic screening is provided on the momentum 

scale g'^T. 

After all, in order to obtain an estimate of the transverse interaction rate we consider the nearly static limit (cu <C q) 
of (0) l25|: 
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2q(q2 + 3r7i2)2 ' 



Pt{oJ < g) = 



4(j6 + (37rm2w/2)2- 

(6) 

Inserting into the integrations can be done exactly, leading to 

- Cf as T, 
T^^CFasT \n-^ 

Ois 

(7) 



We observe that the nearly static limit is a good approximation (within 10%) for the longitudinal rate. The logarithmic 
term of the transverse rate comes from assuming a hypothetical infrared cutoff of the order g^T . For example, the 
magnetic screening mass or the interaction rate itself i.e., the imaginary part of the quark propagator p2| , p3| , 
have been suggested as such an infrared regulator. The coefficient k cannot be calculated using the Braaten-Pisarski 
technique but must await the developement of non-perturbative methods at finite temperature for dynamical quantities 

13- 

For a rough estimate we propose 

r, = r;, + r;~ (2±i)C;^a,T, (8) 

assuming as — 0.3 under the logarithm and k^2. A justification for the latter assumption may come from choosing 
fiD as an upper limit and rrimag as a lower limit for the integration over q. The gluon interaction rate, corresponding 
to the scattering processes shown in Fig.lc and Id, is given by Tg = {CA/CF)Yq — (9/4) F^. Because it is also 
independent of the number of flavors, it holds for the QGP as well as for the pure gluon gas. 

Compared to the real part of the dispersion relation for thermal partons, w ~ fc ~ 3T, the damping rate 7 = F/2 
defined by the imaginary part of the self energy, is not really small, 7^ ~ iasT for realistic values of as- This 
anomalously large damping |^,^ indicates that interactions in the QGP are important, at least for temperatures 
not too far above the phase transition in accordance with the recently emerging picture of the QGP ]42| . 

An alternative method to the Braaten-Pisarski method based on (0) is given by inserting the t-channel diagrams 
of Fig.l in the —t = —{P — P'Y <^ s = (P + K)^ approximation into (P, where the effective gluon propagator A; t is 
used instead of the bare one ji^ . This method has been shown to be equivalent to the self energy calculation (Eq. (H) 
to (^) in the case of the energy loss of a massive fermion in a hot plasma 

Recently, Pisarski proposed an empirical way of including the magnetic mass and an imaginary part of the fermion 
self energy in the parton damping rates |Q . In the case of a hard, massless quark, using the magnetic mass of [^,0 , 
F^ ~ 0.8 Cf as T follows from his investigation assuming as — 0.3 under the logarithm of the transverse part. This 
result is smaller than the estimate (^), since the transverse part of it turns out to be negative for realistic values of 
the coupling constant. 

Next we discuss a much simpler, widely used approximation (see, for example, Ref. [H-p|) based on the naive 
perturbation theory i.e., bare propagators and vertices, where the Debye mass — 3mg is simply introduced by 
hand as an infrared regulator into the gluon propagator. It should be noted that in this case the Debye mass also 
cuts off the magnetic divergence without justification. Then the interaction rate is easily calculated from 

d^k s f , da _ , 

where the parton momentum densities are given by Pq{k) = 12 Nf nF{k) for quarks plus antiquarks and Pg{k) = 
16nB{k) for gluons, respectively. The cross sections in the small momentum transfer limit (— t <^ s) containing the 
Debye mass read 



da 2'na 



C u , ..2,2 ^ (10) 



where the color factor — 4/9 for quark-quark scattering, C = 1 for quark-gluon scattering, and C = 9/4 for 
gluon-gluon scattering. In contrast to the complete calculation (to leading order in the coupling constant) using the 
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Braaten-Pisarski method, the result depends weakly on the number of flavors. In the case of two flavors {Nf = 2) we 
find 



r,(iV/ = 2)~l.la, T, 
Tg{Nf ^2)^2.5asT 



(11) 



and in the pure gluonic case 

r3(A^/ = 0)~2.2a,r. (12) 

Comparing with (||) the simple approach leading to (|l^) or (|l2h appears to underestimate the interaction rates by 
about a factor of two. However, it agrees with Pisarski's result [[4^. 

Finally, we discuss the consequences of the present results for typical values, T — 300 MeV and as = 0.3, expected 
at RHIC and LHC. Using (H) we arrive at relaxation times 

Tg ~ (0.5 ±0.3) /m/c, 
Tq ~ (1.0 ±0.5) fm/c 

(13) 

indicating a rapid thermalization of the QGP and the maintenance of the local thermal equilibrium during the 
expansion phase by comparing with a typical expansion time T^xpan > 1 fm/c Q. Similar results have been found 
from Monte Carlo simulations of ultrarelativistic heavy ion collisions Further more the results ( |l^ ) support 

the prediction of a two-stage equilibration i.e., there is first a thermal equilibrium of the gluonic component before a 
complete thermalization is achieved later on because of the stronger interaction of the gluons compared to the quarks 



III. TRANSPORT INTERACTION RATES 



The relevant physical quantities, such as mean free path and equilibration time, in a plasma with long range 
interactions are described rather by the transport interaction rates than by the ordinary ones [23[|3C|-|3^ . The transport 
rates are obtained from the latter by introducing a weight 1 — cos 9 under the integrals of (|^) or (Q), defining T trans = 
/ dT (1 — cos 9). Here 9 denotes the scattering angle in the center of mass system: cos 9 — (p-p')/(pp') = l + 2t/s. Thus 
the transport factor 1 — cos9 = —2t/s = 2q^ (1 — uP' jq^^j s is proportional to the square of the momentum transfer 

. This additional factor changes the infrared and ultraviolet behavior of the interaction rate completely. The 
transport interaction rate behaves like V trans ^ J dq/q m naive perturbation theory. Therefore the soft as well as the 
hard momentum transfer regimes contribute to T trans ■ This is very similar to the energy loss of a charged particle in a 
relativistic plasma, where an additional factor uj^ appears compared to the interaction rate |^J^ . Quantities which are 
logarithmically infrared divergent in naive perturbation theory turn out to be finite using the Braaten-Pisarski method 
for soft momentum transfers (dynamical screening). Such quantities can be calculated using the method proposed 
by Braaten and Yuan |Q. According to this, introducing a separation scale q*, the soft and hard contributions are 
calculated separately. For the soft contribution (q < q*) resummed propagators and vertices have to be used, whereas 
bare Green's functions are sufficient for the hard one {q > q*). Assuming gT <C <C T the otherwise arbitrary scale 
q* drops out at the end by adding the soft and the hard contribution, reflecting the completeness of the effective 
perturbation theory. In the following the Braaten- Yuan method will be used for computing the transport interaction 
rate following the example of the energy loss |^ as close as possible. 

We will present the calculation of the gluon transport rate in a pure gluon plasma {Nf = 0) in detail, quoting only 
the results for the quark and gluon transport rates in a QGP of two active flavors afterwards. The soft contribution 
follows from (|^) introducing the transport weight 1 — cos 9 under the integral and using q* as an upper limit for the 
q-integration: 



pi{uj,q) + ( 1 - ^ ) Pt(^,g) 



(14) 



The integral over q can be done exactly, while the w-integral has to be evaluated numerically. This has been done 
already in Ref. (l3) yielding {ml = Ana^T'^/S) 



5 



g,trans — g ^9 





- 1.379 







ln(-^ 1-2.811 



(15) 



In Ref. [|3| the hard contribution has not been computed. However, from general arguments (g*-cancellation) we 
know that the hard contribution has to be of the form rg°^^„j = B [ln(r^/g*^) + Ahard], where B = 247ra^r^/s and 
the constant Ahard has to be determined from a detailed calculation of the hard contribution. Thus, if we are only 
interested in a logarithmic accuracy (leading logarithm approximation) valid in the weak coupling limit, we end up 

with m 

Tg.trans = In . (16) 

s as 

However, since the strong coupling constant is not small for realistic values expected in ultrarelativistic heavy ion 
collisions, the knowledge of the coefficient under the logarithm is essential. Therefore the constant Ahard of the hard 
contribution has to be determined. Before we turn to this, let us note that the ln(l/ck5)-term in dlfl) arises from the 
sensitivity of Ttrans to the momentum scale gT, while the logarithmic term in the interaction rate(R) comes from a 
sensitivity to the scale g^T |^^. Thus the transport interaction rate can be calculated completely to leading order in 
the coupling constant using the Braaten-Pisarski method in contrast to F. These entirely different properties of the 
two different kinds of rates originate, of course, from the additional factor in Ttrans- 

For the hard contribution {q > q*) it is sufficient to use the bare gluon propagator. However, in contrast to the soft 
contribution the —t<^s approximation does not hold any more and the u- and s-channel diagrams of Fig.l cannot 
be neglected. The corresponding amplitudes have been calculated a long time ago in Feynman gauge Since 
the on-shell matrix elements are gauge invariant, regardless of the momenta integrated over, there is no problem in 
adopting those results although the soft contribution has been evaluated in Coulomb gauge Pj5l| . 

However, a new difficulty arises now from the divergence of the w-channel contribution for u = [P — K')"^ —>■ i.e., 
—t s. (For massless particles s + t + u = holds.) The u-channel divergence can be regulated in the same way as the 
t-channel singularity, if we choose the transport factor (sin6')^/2 = 2tu/s'^ instead of 1 — cos0. This choice is justified 
because the transport weight 1 — cos 9 has been introduced only for small scattering angles pO| ] for which it is identical 
to (sin 0)^/2. The latter also restores the t-u-channel symmetry in the transport cross sections of quark-quark and 
gluon-gluon scattering. Furthermore, the shear viscosity coefficient beyond the relaxation time approximation also 
contains a factor sin^ 6 ^3 16 3^,Q. Finally, parton collisions with scattering angles near 0° as well as 180° are less 



important for achieving an isotropic momentum distribution (thermal equilibrium) indicating the physical significance 
of a transport rate defined by a weight proportional to sin^ 9 instead of 1 — cos 9 for the equilibration process. 

For calculating the hard contribution of the transport interaction rate in a gluon gas we start from ( l|) modified to 
gluon-gluon scattering and including the factor (sin(?)^/2: 

1 f cPv' f (Pk 

r 9.rarM ^ 0) ^ - J ^ [1 + n^p')] / n,(fc) 

-[1 + UBik')] i2nf 5\P + K-P'-K') 16 {\M{gg -> gg)\^) 



(27r)32fc' 



(17) 



where the matrix element contains the scattering diagrams of Fig. Id. While the hard contribution of the energy loss 
of a heavy fermion with mass M » T could be calculated exactly , this is not possible for (^ . Therefore we 
propose the following approximations: First we assume UBip') — nsip) and nsik') ~ nsik). These simplifications 
hold as long as g = |p — p'l = |k' — k| is not too large or —t is not of the order of s. This assumption may be justified 
because the transport factor {sin 9)"^/ 2 cuts off those momenta effectively. While we will neglect nsip) because of 
(p) ~ 3T, we do not set nsik) = 0. As a matter of fact, the Bose enhancement factor 1 -I- nB{k) is important for the 
exact matching of the soft and hard parts i.e., for the cancellation of q*. 

Using the definition of the differential cross section ||5|] ( [l7| ) may now be written as 



T,,trar^s{Nf = 0) = / ^ 16nB{k) [1 + ^^(fc)] J dt ^. (18) 
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The fc-integration over the Bose distribution functions gives a factor 8T^/3, compared to 16({3)T^ /tt'^ = 1.95 
neglecting the Bose enhancement factor. The differential cross section for gg — > gg scattering according to the 
diagrams of Fig. Id are taken from Ref. M] 



9g4 



us 
1^ 



st 

1 1 2 



tu 



where a factor J/2 has been included to account for the identical particles in the final state 
follows from ( [l8|) by restricting the ^-integration from — s to —q*^. Assuming s ^ q*^ , ( |l8 
in 



247ra?r3 



In 



T^J 15 



(19) 

^ . The hard contribution 
together with (uM results 



(20) 



In order to proceed with the calculation we replace the Mandelstam variable s under the logarithm by its average 
thermal value (s) = 2{p) (k) — 14.59 for gluon momenta (p) — (fc) — 2.701 T. Then we arrive at 



^g%insi^f - 0) - 



247ra2T3 



1.414 



(21) 



Adding up the soft and hard contributions, (p^) and (p^), the separation scale q* drops out as required: 



Tg, trans (A^/ — 0) 



247ra?T3 



In 



1.397 



247ra2T3 , 0.25 
In 



5.2 a^T In 



as 
0.25 



(22) 



where we have used s ~ (s) in the last equation. 

The result ( p^ ) may be compared to the simple approximation of using bare propagators with the Debye mass 
as infrared regulator. Since the degree of the infrared singularity is only logarithmic this amounts approximately to 
using (|l8|), where the separation scale q*^ is replaced by = inagT^ in the upper limit of the t-integration. In this 
way we find 



r 



g.trans 



In 



0.33 



(23) 

may be obtained by using an 



Comparing to (p2[), one realizes that the - to leading order in as exact - result 
effective infrared cut-off of 1.32 instead of 

In the case of a QGP with two flavors we have to consider the diagrams of Fig.la-c in addition. Modifying (|l^ ) 
to these processes, the corresponding calculations are a little bit more involved than in the purely gluonic case. For 
example, we have to be careful about the flavors of the final state e.g., the flavors of the final state quarks of the 
u- and s-channel diagrams of Fig. la and lb have to be identical. Furthermore there is no i-u-channel symmetry in 
the quark-gluon scattering process. After all the u-channel singularity is cancelled by the transport factor (sin 0)^/2 
because it is only logarithmic in this case, rendering the use of an effective quark propagator unnecessary. Also the 
gluon "mass" mg depends now on Nf and, finally, in the case of quarks (or antiquarks) with momenta K and K' we 
have to replace in (|l8|) the factor 16 by 6A^/, nsik) by np{k), and 1 -f ^^(fc) by 1 — npik). The soft contribution for 
the quark rate follows from (|lj) by replacing Ca by Cp- Putting everything together we end up with 

, / Nf\ 247ra,T3 , 0.17 
2=1 + -^ — In 



g,tra7is 



- q,trans 



(Nf 



(Nf 



6 



6.6 ag T In 
6 

2.5 a? T In 



0.17 

as 

'S2TTasT^ 
0.14 



In 



0.14 



(24) 
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where s' = {\ + N f /&) / {1/ Sgg + N f /Qsgq) and s" = {\ + N f /&) / {1/ Sgq + Nf /Qsqq). Here Sgg is the Mandelstam variable 
for two gluons, Sgq for one gluon and one quark, and Sqq for two quarks in the initial state, leading to (s') = 15.13 
and (s") = 17.65 

Applying the transport interaction rates ( p2|) and (^) to ultrarelativistic heavy ion collisions we encounter a serious 
problem: The results (^2|) and (|2^), obtained to lowest order perturbation theory, become negative for realistic values 
of as > 0.2. Hence we cannot draw any conclusions regarding thermalization times at RHIC and LHC from (|24|). 
However, the validity of the Landau collision integral containing the transport cross section beyond the logarithmic 
approximation depends on the condition that the characteristic length / over which the distribution function varies 
significantly must be large compared with the screening length 1/^_d i.e., I fio ^ ^ |0, which is not fulfilled for 
l/l ~ T and /i^i ~ Ay/a^T. Thus the physical significance of the transport interaction rate is somewhat obscure in 
the QGP. 

The unphysical, negative results for the transport rates for realistic values of ctg arise from the separate calculation of 
the soft and hard contributions, which works only in the weak coupling limit g <^ I. In the soft and hard contributions 
the assumption T ^ q* ^ gT is essential for achieving the cancellation of q* . Of course, this assumption cannot 
be fulfilled for g~l rendering (^5|) and ( pT| ) negative, although each contribution is positive by itself without any 
restriction on g*, since they are equivalent to integrals over squares of magnitudes of matrix elements. A similar 
problem occured in ( p3|) by using /i^j as an upper limit for the ^-integration instead of a regulator in the gluon 

propagator and assuming s 3> ^u, which does not hold for g^\ any more. 

The problem of an unphysical, negative result already appeared in the computation of the energy loss of a heavy 
quark in the QGP There it was argued that by using the effective perturbation theory for the entire momentum 
range this problem may be circumvented, increasing, however, the complexity of the calculation drastically. It is not 
possible simply to use the effective gluon propagator for the whole integration range, since it is well defined only for 
tj, g <C T. Gauge invariance and completeness in the order of the coupling constant then demand the use of effective 
vertices and quark propagators in addition, corresponding to including higher orders in g H,^. Thus, in order to 

guarantee positive results and gauge invariance at the same time for values of the coupling constant (7~1, one has to 
go beyond the lowest order in g. In the weak coupling limit, on the other hand, the lowest order contribution, here 
ln(l/as), is sufficient. Furthermore such a unphysical behavior for realistic values of the coupling constant has 
also been observed in the calculation of the gluon plasma frequency beyond leading order |Q. In all of these cases, 
a negative result shows up, if 5 exceeds a critical value of about 1. 

Finally, we will discuss the implications of the transport interaction rates for the viscosity of the QGP 

[^,|6 32 4^ ,[5^-^. In Ref. [^|6| the shear viscosity coefficient -q has been calculated from 



15 rr 



(25) 



where i = q,g, is the energy density of the quarks and gluons in the QGP, and — 2Ti^trans jlSl . (The factor of 
two comes from using the weight sin^ 6 instead of (sin 0)^/2 in the definition of the shear viscosity coefficient p3[.) 
Inserting (EJ) into (psj) we find 



rp3 



0.11 



0.37 



ln(0.17/as) HO.U/as) 



(26) 



This should be compared to the most elaborate calculation of ij (albeit in the leading logarithm approximation) using 
a variational method for the Boltzmann equation resulting in 7] = 1.16r^/[as ln(l/Q!s)] |4^. For as < 0.1 both results 
are comparable. 



IV. CONCLUSIONS 



Before presenting our conclusions, we discuss the validity of the approximations used here. The distinction between 
soft and hard momenta and the cancellation of the separation scale q* , which are the foundations of the Braaten- 
Pisarski and the Braaten-Yuan method p8| , p8[ , rely on the weak coupling limit assumption g 1. In contrast, realistic 
values of as > 0.2 imply g > 1.5. Since as is expected to decrease only logarithmically with increasing temperature, 
g ^ 1 is not even fulfilled at the extreme temperature of the Planck scale. On the other hand, the Braaten-Pisarski 
method is nothing but an improvement of the usual perturbation theory at finite temperature which should work at 
temperatures above twice the critical according to comparisions with lattice QCD [^ , ^ . Furthermore, comparing 
the effective gluon "mass" m^, calculated using the non-perturbative Hartree approximation, with perturbative results 
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shows a difference between the both resuhs by less than 30% at g = 1.5 Also in the case of the energy loss of an 
energetic quark in the QGP it has been shown that the result depends only weakly on the assumption gT <C g* ^ T 
. Those observations indicate that the assumption .g <C 1 should be merely regarded as a mathematical trick and 
not as a physical restriction. Therefore we believe that the Braaten-Pisarski method not only provides a consistent 
treatment of QCD at high temperatures taking into account at the same time important physics as collective effects of 
the non- ideal relativistic plasma e.g., screening, but also gives results for realistic situations which are correct within 
about a factor of 2 as long as logarithmic factors hi{const / as) are not too close to zero or negative as it is the 
case for the transport interaction rate. 

We have estimated the ordinary interaction rate of thermal quarks and gluons by using the effective perturbation 
theory of Braaten and Pisarski, for which the use of an effective gluon propagator is sufficient in the case of thermal 
partons. Due to the missing static magnetic screening in the transverse part of the effective gluon propagator and the 
absence of an imaginary part of the quark propagator we still encounter a logarithmic infrared singularity. Assuming 
a reasonable cut-off, a rough estimate has been obtained, Tg = (6.0 ± 3.0) Us T for gluons and — (2.7 ± 1.3) T 
for quarks, which corresponds to relaxation times of the order r = (0.5 ± 0.3) fm/c for gluons and r — (1.0 ± 0.5) 
fm/c for quarks. This indicates a rapid thcrmalization of the gluon component (two-stage equilibration 0|) and a 
maintenance of the local thermal equilibrium during the expansion phase of the possibly formed QGP at RHIC and 
LHC in accordance with computer simulations of ultrarelativistic heavy ion collisions P,^p^] . 

On the other hand, in a plasma with long range interactions as in QCD the physically relevant quantity should 
be the transport rather than the ordinary interaction rate. The transport rate follows from the ordinary one by 
introducing a transport weight containing the scattering angle in the center of mass system. Due to this factor the 
infrared behavior of the interaction rate is completely changed. The transport rate turns out to be finite using the 
Braaten-Pisarski method because dynamical screening suffices now. We have calculated the transport interaction rate 
for thermal quarks and gluons beyond the leading logarithm approximation by decomposing it into soft and hard 
parts according to the prescription of Braaten and Yuan . The soft contribution has been computed by using the 
effective gluon propagator of the Braaten-Pisarski method, while the hard contribution has been treated using bare 
propagators and vertices. 

Compared to the ordinary interaction rate the transport rate is reduced by a factor of Ug caused by the improved 
infrared behavior due to the transport weight. For a QGP of two active flavors Vg.trans — 6.6 a^T ln(0.17/as) for 
gluons and Tq trans — 2.5 T ln(0.14/as) for quarks have been found. The surprisingly small values of the coefficients 

under the logarithm show that T trans is only meaningful for a .,^0.1. Improving the calculation by using the Braaten- 
Pisarski method over the entire momentum range increases the complexity of the calculation enormously, including 
higher orders of g. Hence no statement about the consequences (e.g., thcrmalization times, mean free paths) for 
realistic values of the coupling constant can be given here. For this purpose, at least a calculation beyond the lowest 
order perturbation theory is required. However, the transport interaction rate obtained here suggests that it may 
be much smaller than the ordinary rate. Thus the realization of a local thermal equilibrium in relativistic heavy ion 
collisions seems to be questionable assuming the transport rates to be responsible for thcrmalization, in contrast to 
computer simulations. However, neither in HIJING [|| nor in the parton cascade 1^ transport cross sections for the 
fundamental parton interactions are used. 

Furthermore the shear viscosity, which is proportional to the inverse of the transport interaction rate , has been 
obtained for the first time beyond the leading logarithm approximation. For values of ofg < 0.1, for which the result is 
well defined, it is large and comparable to the ones obtained by using the leading logarithm approximation ||l^,|l^,^ . 
This supports the idea that dissipation cannot be neglected in hydrodynamic descriptions of the expansion phase of 
the QGP in ultrarelativistic heavy ion collisions [|6|. 

Finally, we have compared the rates obtained from the Braaten-Pisarski method, which arc complete to leading 
order in the coupling constant, with the widely used approach of using bare propagators including the Debye mass as 
an infrared regulator. While the latter approximation works well for transport rates and energy losses ||], it seems 
to underestimate the ordinary rates. This observation suggests that the use of the Debye regulator is justified for 
quantities which are logarithmically infrared divergent in naive perturbation theory as the transport rates or the 
energy loss, but might be questionable for quadratically infrared divergent quantities as the ordinary interaction rate. 



ACKNOWLEDGMENTS 



The author would like to thank T.S. Biro for valuable discussions. The work was supported by the BMFT and GSI 
Darmstadt. 



9 



[1] E. Shuryak, Phys. Rev. Lett. 68, 3270 (1992). 

[2] X. N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 (1991). 

[3] K. Geiger and B. Miiller, NucL Phys. B369, 600 (1992). 

[4] T. S. Biro, E. van Doom, B. Miiller, M. H. Thoma, and X. N. Wang, Duke Univ. Report No. DUKE-TH-93-46, 1993 (to 
be published in Phys. Rev. C). 

[5] J. D. Bjorkcn, Fcrmilab Report No. PUB-82/59-THY, 1982 (unpublished). 

[6] B. Svetitsky, Phys. Rev. D 37, 2484 (1988). 

[7] M. H. Thoma and M. Gyulassy, Nucl. Phys. B351, 491 (1991). 

[8] E. Braaten and M. H. Thoma, Phys. Rev. D 44, 1298 (1991). 

[9] E. Braaten and M. H. Thoma, Phys. Rev. D 44, R2625 (1991). 
[10] S. Mrowczyriski, Phys. Lett. B 269, 383 (1991). 
[11] M. H. Thoma, Phys. Lett. B 273, 128 (1991). 
[12] Y. Koike and T. Matsui, Phys. Rev. D 45, 3237 (1992). 
[13] M. H. Thoma, Phys. Lett. B 269, 144 (1991). 
[14] M. Gyulassy and M. Pliimor, Phys. Lett. B 243, 432 (1990). 

[15] M. Gyulassy, M. Pliimer, M. H. Thoma, and X. N. Wang, Nucl. Phys. A538, 37c (1992). 
[16] P. Danielewicz and M. Gyulassy, Phys. Rev. D 31, 53 (1985). 

[17] See, for example, U. Heinz, K. Kajantie, and T. Toimela, Ann. Phys. (N.Y.) 176, 218 (1987). 

[18] E. Braaten and R. D. Pisarski, Nucl. Phys. B337, 569 (1990). 
[19] R. Baier, G. Kunstatter, and D. Schiff, Phys. Rev. D 45, 4381 (1992). 
[20] H. Nakkagawa, A. Niegawa, and B. Pire, Phys. Lett. B 294, 396 (1992). 
[21] E. Braaten and R. D. Pisarski, Phys. Rev. D 46, 1829 (1992). 

[22] T. Altherr, E. Petitgirad, and T. del Rio Gaztelurrutia, Phys. Rev. D 47, 703 (1993). 

[23] V. V. Lebedev and A. V. Smilga, Ann. Phys. (N.Y.) 202, 229 (1990). 

[24] V. V. Lebedev and A. V. Smilga, Phys. Lett. B 253, 231 (1991). 

[25] R. D. Pisarski, Phys. Rev. Lett. 63, 1129 (1989). 

[26] E. Braaten, Nucl. Phys. B (Proc. Suppl.) 23B, 351 (1991). 

[27] C. P. Burgess and A. L. Marini, Phys. Rev. D 45, 17 (1992). 

[28] A. Rebhan, Phys. Rev. D 46, 482 (1992). 

[29] H. Heiselberg and C. J. Pethick, Phys. Rev. D 47, 769 (1993). 

[30] E. M. Lifshitz and Pitaevskii, Physical Kinetics (Pergamon, New York, 1981). 

[31] H. Heiselberg, G. Baym, C. J. Pethick, and J. Popp, Nucl. Phys. A544, 569c (1992). 

[32] S. Mrowczyhski, in Quark-Gluon Plasma, edited by R. Hwa (World Scientific, Singapore, 1990). 

[33] F. Reif, Fundamentals of Statistical and Thermal Physics (McGraw-Hill, New York, 1965). 

[34] C. J. Pethick, G. Baym, and H. Monien, Nucl. Phys. A498, 313c (1989). 

[35] H. A. Weldon, Phys. Rev. D 28, 2007 (1983). 

[36] V. V. Klimov, Zh. Eksp. Teor. Fiz. 82, 336 (1982) [Sov. Phys. JETP 55, 199 (1982)]. 
[37] H. A. Weldon, Phys. Rev. D 26, 1394 (1982). 
[38] R. D. Pisarski, Physica A158, 146 (1989). 

[39] T. A. de Grand and D. Toussaint, Phys. Rev. D 25, 526 (1989). 

[40] T. S. Biro and B. Miiller, Duke Univ. Report No. DUKE-TH-92-42, 1992 (to be pubhshed in Nucl. Phys. A). 

[41] M. H. Thoma, Mod. Phys. Lett. A7, 153 (1992). 

[42] See, for example, T. S. Biro, Int. Jour. Mod. Phys. 1, 39 (1992). 

[43] G. Baym, H. Monien, C. J. Pethick, and D. G. Ravenhall, Phys. Rev. Lett. 64, 1867 (1990). 

[44] R. D. Pisarski, Brookhaven Report No. BNL-P-1/92, 1992 (unpublished). 

[45] K. Geiger, Phys. Rev. D 46, 4965, 4986 (1992). 

[46] B. Miiller and X. N. Wang, Phys. Rev. Lett. 68, 39 (1992). 

[47] T. S. Biro, B. Miiller, and X. N. Wang, Phys. Lett. B 283, 171 (1992). 

[48] E. Braaten and T. C. Yuan, Phys. Rev. Lett. 66, 2183 (1991). 

[49] B. L. Combridge, J. Kripfganz, and J. Ranft, Phys. Lett. 70B, 234 (1977). 

[50] R. Cutler and D. Sivers, Phys. Rev. D 17, 196 (1978). 

[51] H. A. Weldon (private communication). 

[52] J. D. Bjorken and S. D. Drell, Relativistic Quantum Mechanics (McGraw-Hill, New York, 1964). 

[53] H. Schulz, Hannover Univ. Report No. ITP-IJH 8/93, 1993 (unpublished). 

[54] A. Hosoya and K. Kajantie, Nucl. Phys. B250, 666 (1985). 

[55] R. Horseley and W. Schoenmaker, Nucl. Phys. B280, 768 (1987). 

[56] S. V. Ilyin, O. A. Mogilevski, S. S. Smolyansky, and G. M. Zinovjev, Phys. Lett. B 296, 385 (1992). 
[57] See, for example, B. Petersen, Nucl. Phys. A525, 237c (1991). 
[58] M. Gao, Phys. Rev. D 41, 626 (1990). 

[59] J. Kapusta, P. Lichard, and D. Seibert, Phys. Rev. D 44, 2774 (1991). 



10 



FIG. 1. Lowest order Feynman diagrams for qq — » qq (a), qq — > qq (b), qg — » qg (c), and gg — » gg (d) scattering. 

FIG. 2. Lowest order quark self energy contributions to the quark interaction rate using naive perturbation theory (a) and 
using the Braaten-Pisarski method (b). 
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